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Abstract 

Available experimental data on the exclusive pd — > pnp reaction at 585 MeV show 
a narrow peak in the proton-neutron final-state interaction region. It was supposed 
previously, on the basis of a phenomenological analysis of the shape of this peak, 
that the final spin-singlet pn state provided about one third of the observed cross 
section. By comparing the absolute value of the measured cross section with that 
of pd elastic scattering using the Faldt- Wilkin extrapolation theorem, it is shown 
here that the pd — > pnp data can be explained mainly by the spin-triplet final state 
with a singlet admixture of a few percent. The smallness of the singlet contribution 
is compatible with existing pN — > piVvr data and the one-pion exchange mechanism 
of the pd pnp reaction. 
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Recently, the NN — > NNn reactions with the formation of a spin-singlet NN 
pair in the final state have received a renewed interest. Analyzes of the experi- 
mental data obtained at COSY [1], CESLIUS [2] and LAMPF [3], employ- 
ing the largely model- independent approach of Ref . [4] , show that the singlet 
channel is strongly suppressed in the pp pmi^ reaction at proton kinetic 
energies between 300 and 800 MeV [5-7]. Direct measurements of the singlet 
channel in the reaction pp ppn'^ at RCNP [8] and CELSIUS [9] at 300-400 
MeV indicate a singlet-to-triplet (s/t) ratio of about 1% in coUinear kine- 
matics, which increases up to ~ 10% as the cm scattering angle approaches 
90°. The dominance of the triplet state can be related to the excitation of a 
A-isobar in the intermediate state [7] . 

The measured pion production cross section in pp collision allows one to esti- 
mate qualitatively the s/t ratio in the dcuteron breakup reaction pd — > {pn}p, 
when the quasi-bound {pn} pair is observed in the final state interaction 
(fsi) region and the second proton is detected at large cm scattering angle 
{9* > 90°). It is well known that in backward elastic pd scattering pd dp 
the triangle diagram of one-pion exchange with the subprocess pp — > dir'^ con- 
siderably contributes in the A-region [10]. This mechanism describes well the 
energy dependence of the pd — > dp cross section at 9* = 180° and, in addi- 
tion, explains the qualitative agreement between the proton vector analyzing 
power Ay from pp — > dir'^ and pd dp, observed in the A-region [11]. If 
one assumes that the triangle diagram with one-pion exchange dominates in 
the break-up pd — > {pn}p at large scattering angles, one would expect in this 
reaction a similar s/t ratio of a few percent, as observed in pp — > prnr'^. For 
the A mechanism of the pd pnp reaction, which dominates the one-pion 
exchange triangle diagram, the product of spin and isospin factors yields a. s/t 
ratio of ^ [12]. In contrast, one should expect a higher s/t ratio of about | for 
the one-nucleon exchange mechanism of the deuteron breakup [12]. It was sug- 
gested in Refs. [12-14] to directly measure the singlet channel in the reaction 
pd — > (pp) (0") -1-71(180") with a pp pair of low relative energy Epp — — 5 MeV 
emitted in forward direction and a neutron going backward. Due to a consid- 
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erable suppression of the A-mechanism in this reaction [12] other mechanisms, 
more sensitive to the short-range structure of the deuteron, are expected to 
become important [15]. 

Recent experimental data on the deuteron breakup reaction dp — > pnp with 
two outgoing nucleons in the fsi region were obtained at Saclay [6] at T^ — 1.6 
GeV in semi-inclusive kinematics and at Dubna [16] at = 2-5 GeV. Earlier, 
a kinematically complete exclusive experiment had been performed at Space 
Radiation Effects Laboratory (SREL) in Virginia [17] at a proton beam kinetic 
energy oi Tp — 585 MeV, covering a region of low relative neutron-proton 
energy E^p = — 5 MeV outside of quasi-free pA^-kinematic. A clear peak 
was observed in the five-fold cross section at Enp ~ 0. Using the Migdal- 
Watson approximation [18,19], the authors of Ref. [17] described the shape of 
the fsi peak by assuming a s/t ratio of one third, which corresponds to the 
spin statistical weights of the singlet and triplet states. A smaller s/t ratio of 
about 10% was obtained from the data of Ref. [6]. The difference is possibly 
related to the different cm scattering angles of protons {9* ~ 90° in Ref. [17] 
and e* ~ 180° in Ref. [6] ). 

However, the fitting procedure described in Ref. [17] is rather ambiguous since 
the absolute value of neither the triplet nor the singlet cross section is known 
and was arbitrarily introduced. The s/t ratio can be deduced in principle from 
the data, taking into account only the strong difference in shape of the singlet 
and triplet peaks (see, for example, Ref. [1]). Unfortunately, the low resolution 
in Enp and limited statistics in the peak do not allow one to effectively use 
this procedure for the data of Ref. [17]. In this case the knowledge of the 
absolute value of the triplet (or singlet) cross section is necessary in order 
to determine the s/t ratio. The triplet cross section can be calculated in a 
model-independent way in terms of the large angle proton-deuteron elastic 
scattering. Here we employ the approach described in Refs. [4-7] to determine 
the triplet cross section and on this basis reanalyze the data of Ref. [17]. 

The SREL data are shown in Fig. 1 as a function of the detected proton 
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momentum. At energies Enp of about 1 MeV the cross section is strongly 
influenced by the np fsi. The shape of this peak is well described by the 
Migdal- Watson formulae [18,19], which take into account the nearby poles in 
the fsi triplet (t) and singlet (5) pn— scattering amplitudes 

dasit)^FSIsit){k)K\A,^t)\\ (1) 

Here As(t) is the production matrix element for the singlet (triplet) state, K is 
the kinematical factor, and FSIs{t) is the Goldberger- Watson factor [19]. The 
latter can be written in the form 

'''' = 

where i — s,t. The relative momentum in the pn system at the relative kinetic 
energy E^p — k^/rriN is denoted by k, uin is the nucleon mass. The parameters 
a and (3 are determined by known properties of the on-shell A'"A'"-scattering 
amplitudes at low energies: at = 0.232 fm~^, as = —0.04 fm~^, /3t = 0.91 
fm~^, f3s = 0.79 fm~^ [20]. Important new information on the mechanism of 
pd pnp and off-shell properties of the A^A'^ system is hidden in the matrix 
elements ^s(t), in particular in the ratio 
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C=h^. (3) 



One can find from Eqs. (1) and (3) the following parametrization for the full 
singlet plus triplet cross section [7] 

FSh 

where dot is the triplet cross section. The second term in the brackets of 
Eq. (4) corresponds to the singlet contribution. 



da^^t = (1 + C ^) da„ (4) 



Using the Faldt- Wilkin extrapolation [4], which relates the bound and the 
scattering S-wave functions in the triplet state at short pn distances r < 
1 fm, and by taking into account the short-range character of the interaction 
mechanism, one can find a definite relation between the matrix elements of 
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the pd — > {pn}tp and pd — > dp reactions [4,6]. The triplet differential cross 
section in the laboratory system can then be written as 

d^at{pd ^ pnp) 1 PiPl^ P{k'^) da 



dpi dfli dfl2 167r3 pq ma Ei \pl E„ - ps ■ p„£'2 1 dQ: 
where 



ipd^dp), (5) 



fik') = ^^^^ (6) 

is the Faldt- Wilkin factor [6], da/dO.* is the pd pd cm cross section. In Eq. 
(5) s denotes the squared invariant mass of the pd system, is the deuteron 
mass, po is the beam momentum, Ei and pi {i = l,2,n) are the laboratory 
energy and momentum of the i-th nucleon in the final state. The indices 1 and 
2 refer to the protons and the neutron is referred as n. The proton scattering 
angles in the pd — > {pn}p and pd — > dp processes can be related to each other, 
if the difference between the effective mass of the final {pn} system and that 
of the deuteron is disregarded, as suggested in Ref. [7]. The result presented 
by Eq. (5) should i) be valid at low relative energies Enp, ii) be independent 
of the form of the A^A^-potcntial and details of the large-angle prf-scattering 
mechanism, and iii) it automatically includes the fsi effects in the triplet pn 
system. On the other hand, this method cannot be used for small-angle pd- 
scattering since the A'"A'"-scattering and bound-state wave functions are very 
different at large NN distances, i.e. at low transferred momenta. 

The value of the differential cross section da/dQ,* in Eq. (5) at Tp — 590 MeV 
and 9* = 92.7° amounts to (30.4^°;^[Jf^]) //b/sr [21]. The SREL experiment 
[17] was carried out at almost the same scattering angle {62 — 93.95° for 
Enp = 0). Other available data [22,23] give larger values for the pd dp cross 
section under similar kinematic conditions. Therefore, in order to estimate an 
upper limit for the s/t ratio we use here only the data from Ref. [21]. As one 
can see from Fig. la, the triplet cross section calculated using Eq. (5) (dashed 
line) overshoots the experimental points in the central region around Enp ~ 0, 
but agrees with the data for Epn > 3 MeV. However, a sizable effect arises from 
averaging of the theoretical results over the experimental angular acceptance 
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and resolution of the spectrometer. In order to take these into account, we 
have carried out a five-dimensional integration of the cross section from Eq. 
(5) with Gaussian distributions, where smearing parameters — 2.55° for the 
polar angles of and ap/p = 0.015 for the momentum p were used in accordance 
with Ref. [17]. For the azimutal angles 0i and 02 the averaging was carried 
out in the interval A0 = ±0.4° with a rectangular distribution. 

After smearing we obtain good agreement both in the shape and in absolute 
value between one data set (Fig. la) and a pure triplet contribution of the 
final pn pair with a = 0.7. A small singlet contribution, corresponding 
to C = 0.02, does not contradict the data (x^ = 0.9), whereas larger values 
C = 0.05 (x^ = 1-8) and ( — 0.10 (x^ = 4.6) result in too large a cross section 
in the vicinity of Enp — 0. The other data set (Fig. lb), obtained at a different 
magnetic field setting, shows also dominance of the triplet contribution and 
allows a small singlet fraction: ((, x^)= (0.0, 2.4,), (0.02, 2.1), (0.05, 2.3), (0.10, 
4.0). However, in this case the becomes worse. Under assumption of ^ = |, 
made in Ref. [17], the absolute value of the cross section in the region around 
Enp — results by a factor 2.5-3 too high compared with the data. 

The accuracy of the approximation by Eq. (5) is estimated in Refs. [4-7] and 
[14] to be better than 5% for Enp < 3 MeV. This error arises from variations of 
the bound and scattering NN wave functions at short distances for low Enp- 
The error of the pd — > dp input is ^ 9% [21]. The systematic uncertainties 
in the measured das+t are not given in [17], here we assume them not to 
exceed 10%. Combining all uncertainties given above, the rfc^+t in Eq. (4) is 
uncertain within 15%. If the measured cross section given in Fig. la is scaled by 
factors ranging from 0.85 to 1.15, our x^(C) analysis shows that the resulting 
C's for minimum x^ range from -1-0.035 to —0.030 with the corresponding 
uncertainties ranging from ^.q^H to i'lg 035 ' respectively. This implies that 
C and AC are both of the order of a few percent, and thus are substantially 
smaller than the spin-statistical factor of | assumed in Ref. [17]. 

The matrix element squared |M|^ shown in Fig. 3 was obtained in Ref. [17] by 
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dividing the raw data point by point by a Monte Carlo Enp energy distribution, 
that includes the phase space factor. By this procedure the authors of Ref. [17] 
minimized the effects from averaging over the detector acceptance. In contrast 
to the production matrix element \A\, defined by Eq. (1), the complete matrix 
element \M\ contains the fsi. The authors of Ref. [17] found that the spin- 
statistical fraction of the singlet of | describes the measured data. However, 
the experimental data contain considerable uncertainties. Therefore, according 
to our calculations, they do not constrain the singlet fraction strongly enough. 
As can be seen from Fig. 3, values ( — 0.05 and 0.30 allow one to fit the 
experimental data equally well (x^ = 1.4 and — 0.9, respectively ), if the 
absolute value of the matrix element |Mp, not given in Ref. [17], is treated as 
a free parameter. The small value of (, which we found from the cross section, 
is compatible with the value ( = Q.19~^q^, resulting from our analysis of the 
X^(C) distribution for the |Mp data. 

To improve the sensitivity to the s/t ratio using the extrapolation theorem of 
Ref. [4,6], the ratio of the pd — > pnp and pd dp cross sections has to be 
established better by a measurement of both reactions in the same experiment. 

— * 

A new measurement of the pd — > pnp reaction at the ANKE spectrometer of 
the proton synchrotron COSY-Jiilich will put more stringent limits on the s/t 
ratio by detecting both protons in the forward-forward or forward-backward 
directions at beam energies Tp = 0.5 — 2.5 GeV [15]. 

In conclusion, by comparing the pd pnp cross section at 585 MeV with that 
of pd dp on the basis of scattering theory, we found that the final state 
spin-triplet contribution is dominant allowing a singlet contribution of a few 
percent. This result is in agreement with existing experimental data on the s/t 
ratio in the reaction pN — > pNir and supports the dominance of the triangle 
diagram with the subprocesses pN — > pN tt in the reaction pd pnp. 
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U. and V. K.) gratefully acknowledge financial support and warm hospitality 
at IKP of Forschungszentrum Jiilich. This work was supported in part by the 
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Fig. 1. Experimental cross section (points) of the pd —>■ pnp reaction from Ref. 
[17] at beam energy 585 MeV and proton laboratory scattering angles 6i = 41°, 
82 = 61° as function of the proton momentum in comparison with our calculations, 
a) The pure triplet contribution calculated with corrections taking into account the 
experimental resolution (full line) and without (dashed), as explained in the text. 
The upper scale shows the relative energy (in MeV) of the pn-pair for 61 = 4:1° . b) 
The same observable as in a) but for another magnetic field setting, compared with 
calculations including the corrections for different s/t ratios C = 0.0 (full line), 0.02 
(dashed), and C = 0.05 (dotted). 
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Fig. 2. The squared matrix element, as obtained in Ref. [17], for arbitrary normal- 
ization is well described by C = 0.05, ^ = 1.4 (full line ) and C, = 0.30 = 0.9 
(dashed). 
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